† Corresponding author. E-mail:
Project supported by the Xinjiang Natural Science Foundation of China (Grant No. 2016D01C003).
An exact solution of a single impurity model is hard to derive since it breaks translation invariance symmetry. We present the exact solution of the spin-1/2 transverse Ising chain imbedded by a spin-1 impurity. Using the hole decomposition scheme, we exactly solve the spin-1 impurity in two subspaces which are generated by a conserved hole operator. The impurity enlarges the energy deformation of the ground state above a pure transverse Ising system without impurity. The specific heat coefficient shows a small anomaly at low temperature for finite size. This indicates that the impurity can tune the ground state from a magnetic impurity space to a non-magnetic impurity space, which only exists for spin-1 impurity comparing with spin-1/2 impurity and a pure transverse Ising chain without impurity. These behaviors essentially come from adding impurity freedom, which induces a competition between hole and fermion excitation depending on the coupling strength with its neighbor and the single-ion anisotropy.
The low dimension spin system had revealed many interesting spin orders at low temperature, such as spin-liquids without long range order, spin excitations with energy gap, ferromagnetic order or antiferromagnetic order, and spin glasses.[1–6] These spin orders may demonstrate significant change when a magnetic or non-magnetic impurity is imbedded in the spin chain.[7–11] An antiferromagnetic Heisenberg spin S = 1/2 chain with one impurity is a typical model under frequent theoretical and experimental investigations. Some experimental implementation of these impurity models had been found in the spin-Peierls chain CuGeO3 and spin-ladder SrCu2O3 in which the magnetic atoms Cu are substituted by a few percent of non-magnetic Zn or Mg.[12–14] These non-magnetic impurities result in a collective magnetic order freezing phase at low temperature.[15,16] When a few magnetic Cu atoms in SrCuO2 are substituted by Ni atoms as impurity, the low-energy spin excitation spectrum would be depleted and a spin-pseudogap is open. This character is captured by the one-dimensional Heisenberg model with site defects.[17] Theoretical researches suggests that impurity has a strong influence on electron transportation. Kane and Fisher showed that even a small impurity atom could reflect electrons in a one-channel Luttinger liquid and make the one dimensional link completely insulating.[18]
An exact solution of quantum spin models provides the most accurate predictions of physical phenomena.[19] In some special cases, the exact solution of a Heisenberg chain with an impurity problem was studied by Andrei and Johannesson using Bethe’s Ansatz.[20] They modeled a spin-1/2 Heisenberg chain embedded by impurity atoms with spin S > 1/2. Later on, Lee and Schlottmann generalized this model to impurity with arbitrary spins, and exactly diagonalized the model Hamiltonian to reach its thermodynamics behavior.[21] In the Heisenberg model, all the three components of spin are involved in physics. To simplify the impurity problem in the quantum spin model, a simpler quantum spin model, the quantum transverse Ising model (TIM), also attracted broad interest.[22] To some extent, TIM could be viewed as a special case of the Heisenberg model. The quantum spin Ising model could be exactly solved by a hole decomposition scheme (HDS) method.[23–25] We also used the hole decomposition to derive the exact solution of the spin-1 quantum Ising model.[26,27] The impurity system mentioned above includes many impurity atoms, while the single impurity model is more complicated since it breaks translation symmetry.[28,29] An exact solution of the single impurity model may have potential applications in quantum information for manipulating qubits.
We consider spin-1 impurity embedded in the spin-1/2 Ising spin chain, and use the HDS method and Jordan–Wigner transformation to map the model into the usual spin-1/2 transverse Ising model.[23] The exact solution we obtained reveals many fine structures of quantum critical behavior induced by impurity. This method also can be used to deal with two or three dimensional spin chain with spin-1 impurity.
We consider N spin-1/2 particles coupled to a spin-1 impurity. The neighboring spin-1/2 particles couple to each other by Ising model coupling rule
(1) |
The eigenstates of this model could be expressed in the Hilbert space of
(2) |
For the magnetic impurity space p = 0, the transverse Ising model can be reduced to a quadratic fermion model by Jordan–Wigner transformation[31]
(3) |
(4) |
(5) |
(6) |
In non-magnetic impurity space, the Hamiltonian Eq. (
(7) |
(8) |
This equation can also be obtained from the secular equation Eq. (
The lowest energy level represents the ground state. We compare the lowest energy level in both p = 0 and p = 1 space
(9) |
(10) |
The quantum critical behaviors of the system are related to the impurity excitation, which have two ways: the hole excitation gap Δh and the impurity fermion excitation gap Δx. The excitation gap is determined by the minimal energy for creating a hole or an impurity fermion from the ground states to the excited states. The impurity excitation gap Δimp = min{Δh,Δx} depends on the exchange coupling and single-ion anisotropic of impurity, and the competition of Δx and Δh can induce a switch of the ground state between the two subsystems p = 0 and p = 1. For a positive anisotropic interaction Dz > 0, the magnetic impurity subsystem p = 0 has the lowest energy, so the ground state is confined in the magnetic impurity subspace. If the anisotropic interaction is zero Dz = 0, the system trends to a one-hole subsystem if and only if Jimp = 0. For a negative anisotropic interaction Dz < 0, the non-magnetic impurity subsystem p = 1 bears the lowest energy. Figure
Since the transverse Ising model has an exact solution of ground state energy, we can also exactly find the ground state energy of an impurity imbedded transverse Ising model. Here we are interested in how much energy arises for the ground state when an impurity is introduced into the Ising spin chain. We compute the impurity-induced deformation of the ground state ΔE(p,x) between a pure transverse Ising model and an impurity imbedded transverse Ising model
(11) |
Figure
For the pure transverse Ising model without impurity, the quantum fluctuation due to a small effective magnetic field hx tend to drive the spin chain away from an ordinary magnetic order, while the doping of an impurity into the spin chain enhances the disordered scattering. As shown by the impurity-induced deformation in Fig.
We exactly plotted the competitive relation of the impurity coupling strength and the effective magnetic field at the ground state in Fig.
Specific heat as a thermodynamic quantity can demonstrate the quantum critical point.[33] For example, the specific heat of superfluid 4He has a sharp peak at the temperature T = 2.17 K which marks the quantum phase transition from normal fluid to superfluid. We can exactly compute the specific heat of this impurity model based on HDS. The total partition function of this model is the sum of that of the two sectors, Z = Z(p = 0) + Z(p = 1). As long as we have the partition function, the free energy, entropy, and specific heat can be computed exactly.
Figure
We also calculate the specific heat coefficient at low temperatures for finite size, which is identical to the scaled free energy. From the specific heat coefficient we can observe the thermodynamic anomalous behavior in detail, which should be capable of yielding important information on the nature of the impurity.[34,35] In Fig.
In Fig.
Notice that the impurity above is a spin-1 particle. One may ask if the spin-1 impurity is replaced by a special spin-1/2 particle, whether it would induce a similar anomaly or not. So we calculated S = 1/2 impurity in the same model as the S = 1 impurity. The light blue dash line in Fig.
A closer observation on Fig.
Impurity models are only exactly solvable in some special cases. As the impurity breaks the translational invariant symmetry, it is hard to derive an exact result for a quantum impurity model. We exactly solved the transverse Ising model imbedded with an S = 1 quantum impurity. The Hilbert space of this model is divided into two spaces due to a conserved quantum operator, the hole operator. With this exact spectrum of ground state and excited states, we found that the impurity could enlarge the energy deformation of the ground state above a pure transverse Ising model. The deformation grows approximately linearly with respect to an effective magnetic field. We computed the dependence of the ground state on impurity parameters. There exists a tuned behavior of the ground state, which depends on the effective magnetic field and the impurity coupling strength. We computed the specific heat coefficient to explore the anomalous behavior at low temperature for finite size. The specific heat coefficient shows a sharp peak when the impurity coupling strength becomes weak enough. While this peak is the special transition point induced by spin-1 impurity, there is no such sharp transition point for a pure transverse Ising model and spin-1/2 impurity model. This transition point indicates the switch of the ground state from a magnetic impurity space to a non-magnetic impurity space. Since the magnetic impurity space to a non-magnetic impurity space is determined by the eigenvalue of the conserved quantum operator of spin-1, it is essentially a competition of hole excitation and fermion excitation. As the spin-1 impurity could induce some special thermodynamics, it maybe also has some exotic effect on the quantum entanglement at low temperature. In the future, we shall study the quantum entanglement, geometric phase, and potential topological phases.
[1] | |
[2] | |
[3] | |
[4] | |
[5] | |
[6] | |
[7] | |
[8] | |
[9] | |
[10] | |
[11] | |
[12] | |
[13] | |
[14] | |
[15] | |
[16] | |
[17] | |
[18] | |
[19] | |
[20] | |
[21] | |
[22] | |
[23] | |
[24] | |
[25] | |
[26] | |
[27] | |
[28] | |
[29] | |
[30] | |
[31] | |
[32] | |
[33] | |
[34] | |
[35] |